Lagrangian Statistics and Temporal Intermittency in a Shell Model of Turbulence 
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We study the statistics of single particle Lagrangian velocity in a shell model of turbulence. We 
show that the small scale velocity fluctuations are intermittent, with scaling exponents connected to 
the Eulerian structure function scaling exponents. The observed reduced scaling range is interpreted 
as a manifestation of the intermediate dissipative range, as it disappears in a Gaussian model of 
turbulence. 



In the recent years there has been a great improvement 
in the experimental investigation of turbulence from a La- 
grangian point of view . In the Lagrangian approach 
the flow is described by the (Lagrangian) velocity v{xq, t) 
of a fluid particle initially at position x(Q) = x$. This is 
the natural description for studying transport and mixing 
of neutrally advected substances in turbulent flows. 

One of the simplest statistical quantities one can be in- 
terested in is single particle velocity increments Sv(t) = 
v(t) — v(0) (where, assuming statistical homogeneity, we 
have dropped the dependence on Xq) for which dimen- 
sional analysis in fully developed turbulence predicts 



(1) 



where e is the mean energy dissipation and Co is a nu- 
merical constant. The remarkable coincidence that the 
variance of Sv(t) grows linearly with time is the physi- 
cal basis on which stochastic models of particle disper- 
sion are based. It is important to recall that the "dif- 
fusive" nature of (Q) is purely incidental: it is a direct 
consequence of Kolmogorov scaling in the inertial range 
of turbulence and is not directly related to a diffusive 
process. Let us recall briefly the argument leading to 
the scaling in ([!]). We can think at the velocity v(t) ad- 
vecting the Lagrangian trajectory as the superposition of 
the different velocity contributions coming from the tur- 
bulent eddies (which also move with the same velocity 
of the Lagrangian trajectory). After a time t the compo- 
nents associated to the smaller (and faster) eddies, below 
a certain scale £ are decorrelated and thus at the lead- 
ing order one has Sv(t) ~ Sv(£). Within Kolmogorov 
scaling, the velocity fluctuation at scale £ is given by 
Sv(£) ~ V^/L) 1 / 3, where V represents the typical veloc- 
ity at the largest scale L. The correlation time of Sv(£) 
scales as t{£) ~ t q (£/L) 2 ^ 3 and thus one obtains the scal- 
ing in (0). 

Equation (gj) can be generalized to higher order mo- 
ments with the introduction of a set of temporal scaling 
exponent £(p) 



(2) 



diction = p/2 but one expects corrections to dimen- 
sional scaling in presence of intermittency. 

A generalization of (|l|) which takes into account inter- 
mittency corrections can be easily developed within the 
multifractal model of turbulence HQ. The dimensional 
argument is repeated for the local scaling exponent h, 
giving Sv(t) ~ V(t/To) h ^ 1 ~ h \ Integrating over the h 
distribution one ends with the prediction 



(Sv(t) p ) ~V P dh 



\ph-D(h)+3]/(l-h) 



(3) 



where, in the limit t/ro — > 0, the integral can be esti- 
mated by a steepest descent argument giving finally 



£(p) = min 

h 



ph- D{h) + 3 



(4) 



The dimensional estimation sketched above gives the pre- 



The fractal dimension D(h) is related to the Eulerian 
structure function scaling exponents C(?) by the Legen- 
dre transform [Q ((q) — mmh[qh — D(h) + 3]. The stan- 
dard inequality in the multifractal model (following from 
the exact result ((3) = 1), D(h) <3h + 2 implies for (@) 
that even in presence of intermittency £(2) = 1. This is 
a direct consequence of the fact that energy dissipation 
enters in ([!]) at the first power. Our expression for scaling 
exponents (Q) recovers in a more compact form the pre- 
diction obtained on the basis of an "ergodic hypothesis" 
of the statistics of energy dissipation j|] . 

Recent experimental results [Q have shown that in- 
deed Lagrangian velocity fluctuations are intermittent 
and display anomalous scaling exponents, as predicted by 
the above argument. Despite the relative high Reynolds 
number of the experiments, a true temporal scaling range 
is not observed. Thus the estimation of the scaling ex- 
ponents in (^) can be done only relatively to a reference 
moment (the so-called ESS procedure ||). 

In this work we use a dynamical toy model of tur- 
bulence for investigating scaling (||) and prediction ([|) 
at very high Reynolds numbers. The model is a popu- 
lar shell model of turbulence in which the velocity 
fluctuations of the eddies at the scale £ n = L2~ n are rep- 
resented by a single variable u n (n = 1,2, N). Only 
local interaction among shells are represented and there- 
fore no sweeping effects are present. In this sense shell 
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models are dynamical models of velocity fluctuations in 
a Lagrangian framework, and have been already used for 
the study of turbulent dispersion JTl|] . The equation of 
motion for the complex shell variable u n is 



~~dt 



■ if * $ * 

ikn M„ +2 ti„ +1 - -u n+1 u n _ 1 + 
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-U n -lU n -2 



-vk 2 n u n + f n 



(5) 



where k n = i^ 1 and /„ is a deterministic forcing acting 
on the first two shells only. Shell model (||) is charac- 
terized by a chaotic dynamics with a statistically steady 
state with a constant flux of energy from large to small 
scales. The fluctuations generate by chaoticity induce a 
breaking of the global scaling invariance and corrections 
to the Kolmogorov exponents for the structure functions 
close to the experimental values |T(| . 

Lagrangian velocity in the shell model framework can 
be represented as the superposition of the contributions 
of all the different eddies. Let us thus define 



v(t) = ^2 -Re(u„) 



(6) 



where we have taken, rather arbitrarily, only the real part 
of the shell variables with unit coefficient. From the def- 
inition of the shell model, there is not a precise recipe 
for reconstructing the Lagrangian velocity. More in gen- 
eral, one could think of a representation in which shell 
variables are multiplied by an appropriate wavelet func- 
tions. Of course numerical prefactor, such as Cq in (|l|) 
will depend on the wavelet basis. Nevertheless one ex- 
pects that different choices should not affect Lagrangian 
scaling exponents which are determined by the dy- 
namical properties of the model. 

Previous studies of multi-time correlations in shell 
models of turbulence have demonstrated the existence 
of a set of correlation times compatible with the mul- 
tifractal picture of the turbulent cascade jl3|. This is 
an indication that, as we will see, Lagrangian velocity 
defined as (||) will be affected by intermittency. 

Very long and accurate numerical simulations of the 
shell model (§) with N = 24 shells and 6=1/2 have 
been performed. The energy is injected at a constant 
flux e = 0.01 is the first 2 shells and is removed at the 
smallest shells by viscosity v = 10~ 7 . With these param- 
eters, our simulations correspond to a Reynolds number 
Re ~ 10 8 . For each realization, Lagrangian structure 
functions are computed from the Lagrangian velocity (|^) 
up to the large-scale time To- Average is then taken over 
10 5 independent realizations. 

In Fig. [l] we plot the set of numerically determined Eu- 
lerian structure function scaling exponent together 
with the fractal dimension D{h) reconstructed by means 
of the Legendre transform. We observe strong intermit- 
tency in velocity statistics with scaling exponent which 
deviates from Kolmogorov prediction. The scaling ex- 
ponents are not universal with respect to the particular 



shell model. Model (@) gives a set of exponents which 
are a little more intermittent than, but not far from, 
the experimentally observed exponents |Q: £(2) ~ 0.72, 
C(4) ~ 1.25, C(6) ^ 1.71. We thus expect that the values 
of £ (p) obtained from (Q) using the D(h) of Fig. will be 
directly comparable with real experimental data. 

Figured shows the second-order Lagrangian structure 
function (uh as a function of time. The linear behavior is 
evident (see the inset) even if a long crossover from the 
ballistic scaling at short time (Sv(t) 2 ) ~ t 2 is present. In 
spite of the very high Reynolds numbers achievable in 
the shell model, the extension of the temporal scaling (||) 
is still moderate. For a comparison with the available 
experimental data, in the inset we also plot the result 
obtained from two simulations at lower resolution, with 
Re ~ 2 x 10 6 and Re ~ 10 5 . In the latter case almost no 
scaling range is observable. Despite these limitations, we 
will see that high Re simulations allow the determination 
of the Lagrangian scaling exponents with good accuracy. 

The long crossover in Fig. ^| can be understood in terms 
of intermediate dissipative range as a consequence of the 
fluctuating dissipative scale The smallest time 

at which one can expect scaling ([[]) is the Kolmogorov 
time t v ~ r i?e _(1_ ' l) / (1+ ' l) which fluctuates with h. A 
demonstration of the effects induced by intermittency is 
given by considering a non-intermittent Gaussian model. 

Setting f n = v = 0, (^) becomes a conservative sys- 
tem with two conserved quantitities which depends on 
the value of the S [ jlOj ] . In statistical stationary con- 
dition, the model shows equipartition of the conserved 
quantities among the shell, in agreement with statistical 
mechanics prediction For 5 = l + 2~ 2 / 3 the equipar- 
tition state leads at small scales to Kolmogorov scaling 
(|u„| 2 ) ~ i 2 / 3 with Gaussian statistics. In Fig. || we plot 
the second-order Lagrangian structure function ([!]) for 
the Gaussian model. Both the ballistic and the diffusive 
scaling is clearly observable and the crossover is strongly 
reduced with respect to Fig. ||. 

In Fig. H we plot the probability density functions of 
5v(t) computed at different t in the linear scaling range of 
Fig. H rescaled with their variances. The form of the pdf 
varies continuously from almost Gaussian at large time 
(i ~ t ) to the development of stretched exponential tails 
at short times, similar to what observed in laboratory 
experiments Q. Flatness F grows from Gaussian value 
F = 3 up to F ~ 20 at smallest times. This is an in- 
dication of Lagrangian intermittency, in the sense that 
the Lagrangian statistics cannot be described in term of 
a single scaling exponent. 

In Fig. H we plot the set of Lagrangian scaling expo- 
nents £(p) obtained from a direct fit of temporal struc- 
ture functions. The nonlinear behavior in p confirms the 
presence of Lagrangian intermittency already observed 
from the pdf. We present the result for moments up 
to p = 6 which approximatively corresponds, from (Q), 
to Eulerian structure function of order q = 8. In this 
sense temporal structure functions are more intermittent. 
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Figure |5j shows that the agreement with the multifractal 
prediction (^) is very good up to the moment achiev- 
able with our statistics. What is even more remarkable 
is that our prediction is very close to experimentally de- 
termined exponents. For example we find £(3) ~ 1.31, 
£(4) ~ 1.58, £(5) ~ 1.85, while the experimental data 
give § £ex P (3) = 1.34 ± 0.02, £ ea;p (4) = 1.56 ± 0.06 and 
£exp( 5 ) = 1.8 ±0.2. 

In this work we have investigated the statistical prop- 
erties of single particle Lagrangian velocity in fully devel- 
oped turbulence. A prediction for intermittent scaling ex- 
ponents of Lagrangian structure function is given within 
the multifractal framework. Very high Reynolds number 
simulations in shell model confirm the multifractal pre- 
diction, even if rather small scaling ranges are observed. 
In particular, our simulations show that at experimental 
Reynolds number presently available almost no scaling 
is observable. The reduction of the scaling range in La- 
grangian statistics is interpreted as an effect of the inter- 
mediate dissipative range. A Gaussian, non-intermittent 
version of the shell model confirms this interpretation. 
An important consequence of our findings is that usual 
models of particle dispersion, based on stochastic model 
p6[ , are incorrect in a fundamental sense, and one should 
take into account the modifications due to intermittency. 

This work was supported by MIUR-Cofin2001 contract 
2001023848. We acknowledge the allocation of computer 
resources from INFM "Progetto Calcolo Parallclo" . 
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FIG. 1. Shell model Eulerian structure function scaling ex- 
ponents £(q) determined from over the statistics of 10 5 inde- 
pendent configurations. In the inset we plot the codimension 
3 — D(h) as determined from numerically solving the Legendre 
transform. 
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FIG. 3. Second-order Lagrangian structure function 
(5v(t) 2 ) as a function of time delay t for the equilibrium Gaus- 
sian model. Continuous line is the ballistic behavior t 2 at 
short time. Dashed line represents the linear growth (hi). 




FIG. 2. Second-order Lagrangian structure function 
(Sv(t) 2 ) as a function of time delay t for the simulation at 
Re — 10 s . Continuous line is the ballistic behavior t 2 at short 
time. Dashed line represents the linear growth (|l|). Inset: 
(5v(t) 2 ) compensated with the dimensional prediction st at 



Re = 10 (continuous line), Re 
Re = 10 5 (dotted line). 



2 x 10 (dashed line) and 




8v/<5v 2 > 1/2 

FIG. 4. Probability density functions of velocity dif- 
ferences Sv(t) normalized with the variance at time lags 
t/r = 0.002(D), 0.01(*),0.06(x),0.35(+). Continuous line 
represents a Gaussian. Inset: flatness F = (Sv(t) ) / (8v(t) 2 ) 2 
as function of time and Gaussian value F = 3 (dashed line). 
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FIG. 5. Lagrangian structure function scaling exponents 
£(p) numerically determined by a best fit on (El). The line 
represents the multifractal prediction (^) with D(h) obtained 
4 from Fig. [j]. 



